Profile and scaling of the fractal exponent of percolations in complex networks 
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We propose a novel finite size scaling analysis for percolation transition observed in complex 
networks. While it is known that the cooperative systems in growing network models often undergo 
an infinite order transition with inverted Berezinskii-Kosterlitz-Thouless singularity, the transition 
point is very hard for numerical simulations to find because the system is in a critical state outside 
the ordered phase. We propose a finite size scaling form for the order parameter by using the fractal 
exponent, which enables us to determine the transition points and critical exponents numerically 
for infinite order transitions as well as standard second order transitions. We confirm the validity of 
our scaling hypothesis through the Monte-Carlo simulations for bond percolations in two network 
models; the growing random network model, which exhibits an infinite order transition, and its 
reconstructed network, which shows a conventional second order transition. 
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I. INTRODUCTION 

The science of complex networks has led to new per- 
spectives in the study of the statistical physics [U-IH • Var- 
ious processes, such as percolations, epidemic spreadings, 
spin systems, and coupled oscillators, in complex net- 
works have been actively studied (see [f| and references 
therein). 

An elementary theoretical framework of critical phe- 
nomena in complex networks has been provided in the 
manner of the mean field approximation and local tree 
approximation Particularly, the equilibrium systems 
on the random graphs with arbitrary degree distributions 
P(k) (k denotes degree), so-called configuration model, 
are well described by local tree approximation, where 
critical properties are determined only by the power ex- 
ponent 7 of the degree distribution, P(k) oc fc~ 7 @, @. 
On the other hand, the systems on network models made 
with growth mechanism are not the case. Let us con- 
sider the bond percolation with open bond probability p. 
Several analytical works for the percolations on growing 
network models reveal that the system exhibits unusual 
phase transitions, termed infinite order transition with 
inverted Berezinskii-Kosterlitz-Thouless (BKT) singular- 
ity [(3|-[l8| (see also (l9l - l2]| for the case of spin systems): 
(i) The singularity of the phase transition is infinitely 



weak. When p is larger than the transition point p c , the 
order parameter m = S max (N)/N, where S max (N) is the 
mean size of the largest clusters over percolation trials in 
the system with N nodes, follows 



m oc exp(— a/ (Ap)' 3 ), 



(1) 
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where Ap = p — p c . (ii) Below the transition point, the 
mean number n s of clusters with size s per node obeys 
the power law, 

n s oc s~ T . (2) 

Recent studies [13, HH revealed that this phase is simi- 
lar to the critical phase (also termed intermediate phase 
[22I . I23I ]) observed on (transitive) nonamenable graphs 
(NAGs). NAGs are defined as infinite graphs with pos- 
itive Cheeger constant [HI, Typical examples of 
NAGs arc hyperbolic lattices and trees. The bond per- 
colation on a NAG takes three distinct phases; (i) the 
non-percolating phase (0 < p < p c \) in which only finite 
size clusters exist, (ii) the critical phase (p c i < p < p c -2) 
in which there are infinitely many infinite clusters, and 
(iii) the percolating phase (p C 2 < p < 1) in which the sys- 
tem has a unique infinite cluster. Here an infinite cluster 
means a cluster whose size is of order 0(N a ) (0 < a < 1). 

It has been proven that the phase boundaries of 
amenable graphs including the Euclidean lattices, and 
NAGs satisfy <Pd = Pc2 < 1 and < p c \ < p C 2 < 1, 
respectively [HI, [H[ ■ The scenarios of the phase transi- 
tions of the percolations on regular graphs are summa- 
rized in Table HI If the above-mentioned unusual phase 
on growing networks is regarded as the critical phase on 
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TABLE I: Expected scenario of phase transitions of bond percolations on various graphs. The number of ends of an infinite 
graph G is defined as the supremum of the number of infinite connected components in G\S, where G\S is the graph obtained 
from G by removing arbitrary finite subgraph S and the edges incident to those [2^ . |23| . 
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amenable graph with two ends 


chain 


< Pel = Pc2 = 1 


amenable graph with one end 


d(> 2)-dimensional Euclidean lattice 


< Pel = Pc2 < 1 


NAG with infinitely many ends 


Cayley tree 


< Pel < Pel = 1 


NAG with one end 


enhanced binary tree [24-26] , hyperbolic lattice [2JJ 


< Pel < Pe2 < 1 


tree under growth mechanism 


growing random tree [8, 11. 12. 14. 171 


= Pel < Pc2 = 1 


network under growth mechanism 
hierarchical scale free networks 


growing random network (m-out graph) [11. 28. 291 
Callaway-Hopcroft-Kleinberg-Newman-Strogatz model [6, 7} 
decorated (2,2)-flower [13. 16. 18], Hanoi network [15.1 


= Pel < Pe2 < 1 



NAGs, we may say that the growing networked systems 
add a new scenario; p c \ = and p c2 ^ (last two rows 
in Table |T]). However, it is not known which relation 
(or none of them) holds for complex networked systems, 
and wc should keep the above three phases in mind in 
studying the critical phenomena. This investigation is a 
challenging study because numerical studies so far were 
not conscious of the existence of the critical phase. 



Finite scaling analysis is a powerful method to extract 
the transition point and critical exponents from numer- 
ical data. It is indeed the case for some non-growing 
network models as seen in a finite size scaling analysis 
proposed by Hong et al. (30[. However, to estimate the 
transition point p c2 and critical behaviors of the grow- 
ing networks, we are faced with the following difficulties. 
The transition at p c2 is the one between the percolating 
phase and the critical phase. We should perform a scal- 
ing analysis with data only for p > p c2 because the phase 
below p c2 is critical so that a standard scaling analysis 
cannot be applied [I?}. Second, the transition at p c i is of 
infinite order so that conventional critical exponents do 
not work as normal fitting parameters. 



In this paper, we propose a finite size scaling analy- 
sis to investigate numerically both the growing and non- 
growing networked systems. Our finite size scaling form 
is described in terms of the network size N and the 
fractal exponent, which characterizes the critical phase 
fl7l [HI, H3, . As an application, wc evaluate the frac- 
tal exponent and perform the scaling analysis for the 
percolations both on the random attachment growing 
network (RAGN; a random attachment version of 
the Barabasi- Albert model [311]). which exhibits an in- 
finite order transition with inverted BKT singularities 
[HHHjHilj and its reconstructed network (a random net- 
work having the degree distribution same as the RAGN) , 
which exhibits a second order transition from the non- 
percolating phase to the percolating phase. Our finite 
size scaling analysis works quite well for both cases. 



II. SCALING ARGUMENT 

To investigate the critical phase, it is useful to profile 
the fractal exponent ip 24J. It is defined as 



ip = lim log N S max (N), 

N-too 



(3) 



which mimics df/d for d-dimensional Euclidean lattice 
systems, df being the fractal dimension of the largest 
clusters. In the critical phase, ip is related to r in Eq.© 
as t = l + ip^ 1 , indica ting that ip plays a role of a natural 
cutoff exponent of n s [l7,[l8|. Then, the non-percolating 
phase, the critical phase, and the percolating phase are 
represented with ip = 0, < ip < 1, and ip = 1, respec- 
tively. 

Let us assume a scaling form with ip and N in the or- 
dered phase. First, we suppose the existence of crossover 
size N*(p) at p. When N is much smaller than the 
crossover size, N <C N*(p), the system behaves as if it 
were critical even for p > p c2 ; S max (N) oc iV^' c , while 
S max (N) is proportional to N when N > N*(p). By 
connecting these two limits, we expect a finite size scal- 
ing form for S max (N) as 



t (jyO = ^ a /i( 



N 



N*(p) 



where 



for x 3> 1 



const for x <C 1, 



(4) 



(5) 



ip c being the fractal exponent at p = p c2 . Or cquivalcntly, 

N 

(6) 



S ma *(N)=N*(pf°f 2 



N*{p)J' 



where 



f / s _ \ x for a; » 1 

J2W - S x i, c for X< 1_ 



(7) 



By taking the logarithmic derivative of Eq.([6j) with N, 
we obtain a finite size scaling form for ip as 



iP(N) 



dlnS max {N) 
din AT 



9 



N 



N*(p)r 



(8) 
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FIG. 1: (Color Online) (a) The order parameter m and (b) the 
fractal exponent ip on the RAGN. The numbers of nodes N 
are 2 17 (black-diamond), 2 16 (red-circle), 2 15 (green-square), 
2 14 (blue-triangle), 2 13 (open-diamond), and 2 12 (open-circle). 
The vertical dashed line indicates p C 2 = (1 — l/\/2)/2, and 
the horizontal dashed line indicates ip c = 0.44. 



where 



dlnf 2 {x) 
d\n(x) 



1 for x ^> 1 

ipc for 



(9) 



Similarly, the finite size scaling form for the average size 
S av {N) = X^s m , x s2rLs °f * ne cms ter to which any node 
belongs can be assumed as 



S &V (N) = N^h^ 



N 



N*{p)J' 



where 



h(x) 



x~^ av for x > 1 
const for x <C 1, 



(10) 



(11) 



and ip av = 2ip c — 1, which corresponds with the fractal 
exponent of the mean cluster size at p = p C 2- 

For practical use of Eqs.((4]) and (|SJ), we need to know 
crossover size N*(p). By considering the case N ~ 
N*(p), we relate N*(p) to S ma , x (N) and m(p) as 
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FIG. 2: (a) The order parameter m, (b) the fractal exponent 
^, and (c) the mean cluster size S av on the reconstructed net- 
work. The numbers of nodes N are 2 17 (black-diamond), 2 16 
(red-circle), 2 15 (green-square), 2 14 (blue-triangle), 2 13 (open- 
diamond), and 2 12 (open-circle). The vertical dashed line in- 
dicates p c2 = 2/9, and the horizontal dashed line indicates 
ipc = 2/3. 



In the case of the second order transition, m(p) follows 
m oc (Ap)^ for p > p c2 = p c i. Then the finite size scaling 
form for S max is 

S max (N) = N^hiNiAp)^), (13) 
and that for -0 is 

if) = g(N(Ap) T ^*o). (14) 

This argument also can be applied to the case of the 
infinite order transition. Because m(p) follows Eq.([T]), we 



4 



(a) 

200 
100 
50 

20 
10 
5 




io- 



io- 



10" 



io- 



io- 



10 



(b) 



1. 
0.9 
0.8 
0.7 

0.6 
0.5 



10" 



io- 



10" 7 10" 4 10" 



10 2 



FIG. 3: Finite size scalings for (a) S ma x by Eq. (fT5)) and (b) 



i> by Eq.(|T6)l on the RAGN. N = 2 U (black-diamond), 2 lb 
(red-circle), 2 15 (green-square), 2 14 (blue-triangle), 2 13 (open- 
diamond), and 2 12 (open-circle). The horizontal dashed line 
follows tp c = 0.44. 



obtain finite size scaling forms of infinite order transition 
for S max (N) and tp by substituting Eq.([T2|) into Eqs.(j4]) 
and © as 

S max (N) = N^h{Ne- a '^-^ A ^), (15) 



and 



respectively. 



(16) 



III. NUMERICAL CHECK 

To check the validity of our scaling argument, we per- 
form the Monte-Carlo simulation for the bond percola- 
tions on the RAGN and its reconstructed network. The 
RAGN stochastically generates a graph realization with 
N nodes in the following way. We start at time t = with 
a complete graph with I + 1 nodes. At each time step, 
a new node joins into a graph, and links to I randomly- 
chosen nodes among pre-existing nodes. The network 
growth continues until t = N — I — 1. The station- 
ary degree distribution P(k) of the resulting network is 
P(k) oc (1/(1 + l)) fe . A reconstructed network of the 
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FIG. 4: Finite size scalings for (a) Smax by Eq.|T3]), (b) ip by 
Eq. (|14[) . and Sav by Eq.lJTO]) on the reconstructed network. 
TV = 2 17 (black-diamond), 2 16 (red-circle), 2 15 (green-square), 
2 14 (blue-triangle), 2 13 (open-diamond), and 2 12 (open-circle). 
The horizontal dashed line follows ip c — 2/3. 



RAGN is then obtained by reconnecting its edges at ran- 
dom, keeping the degree of each node. This operation 
generates an uncorrelated network with the same degree 
distributions as the original one. 

The rigorous results of the critical behavior for the 
bond percolation on the RAGN were given by Bollobas 
and Riordan [2^, [2(| ; the transition point p C 2(l), above 
which a giant component having size of O(N) emerges, 
is given by p C 2(l) = (1 — y/l — l/l)/l and for p > p c2 the 
order parameter follows Eq.([T]) with a = tt/2[1(1 — 1)] 1/4 
and/3' = 1/2. The bond percolation on the reconstructed 
network is analyzed by local tree approximation 0, 0). 
The transition point is given by p c = (k) / (k 2 — k) = 2/9. 
Also, it is known that the percolations on uncorrelated 
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networks with P(k) oc k 7 , where 7 > 4, belong to the 
mean field universality class, ft = 1. 

We set I = 2, and N = 2 n ,2 12 ,- • ■ ,2 18 . The number of 
graph realizations is 100000 (10000 for N =2 16 ,2 17 ,2 18 ), 
and the number of percolation trials on each realiza- 
tion is 100. The order parameter m and the frac- 
tal exponent ip on the RAGN arc shown in Fig. [T] 
Here we calculate ip(N) by the difference [In S max (2N) — 
\nS max (N/2))}/[\n(2N) - ln(iV/2)]. As N increases, 
ip(N) converges to a certain value depending on p for 
P < Pc2i while ip(N) for p > p c2 approaches unity very 
slowly with N. In the thermodynamic limit, N — > 00, ip 
grows continuously with p for 0(= p c x) < p < p C 2 up to 
ip c ~ 0.44 at p = p C 2 = (1 — y/l/2)/2, and then jump to 
?/' = l. The same behavior of ip was analytically observed 
on the decorated (2,2)-flower fl8j . 

The fractal exponent is useful not only to confirm the 
existence of the critical phase, but also to determine the 
critical point of ordinary continuous phase transition. In 
Fig. [21 m, ip, and 5 av on the reconstructed networks are 
plotted with respect to p. In this case, ip with several 
sizes cross at (p c , ip c ) ~ (2/9, 2/3). In the limit N — > 00, 
ip = 0(= 1) for p < p c {> p c ), which means the transition 
at p c is the one between the non-percolating phase and 
the percolating phase. The result that a giant component 
at the transition point is of 0(N 2 / 3 ) is also observed on 
the Erdos-Renyi model [32| . 

The finite size scaling results for the RAGN and the 
reconstructed network are shown in Figs. [3] and 01 re- 
spectively. Here fitting parameter ft (and also a for the 
RAGN) is set to the analytically obtained values. We ob- 
serve good data collapses for both cases, confirming the 
validity of our scaling hypothesis. We did not perform 
the finite size scaling of S av for the RAGN. The reason is 
that our numerical result for the RAGN shows ip c < 1/2 
at the transition point p = p C 2, leading to ip av < 0. This 



indicates that the mean cluster size of the RAGN does 
not diverge, and shows a finite jump at p = p C 2- This is 
actually observed on the RAGN [ll| . 



IV. SUMMARY 

To summarize, we have proposed a novel scaling anal- 
ysis for critical phenomena. We have confirmed the va- 
lidity of our scaling via two exactly-solved models; the 
RAGN and its reconstructed network. 

Note that the present scaling form includes the con- 
ventional finite size scaling. For <i-dimensional lattice 
systems, the fractal exponent satisfies ip c = 1 — ft/dv, so 
that Eq.Q reduces to the conventional scaling for Nm 
provided that N/N*(p) = (L/£) d , where L is the lin- 
ear dimension, and v the critical exponent of correlation 
length £. In this sense, N*(p) is regarded as the correla- 
tion volume 

Whether percolations on the other complex networks 
have such a critical phase or not is not known. Also, 
whether other dynamics, e.g., epidemic spreading, cou- 
pled oscillators and so on, have such a critical phase or 
not is an interesting open problem. By answering to these 
questions, the relation between network topology and the 
dynamics thereon will become more clear. We believe 
that the present method will contribute to the further 
analysis. 
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